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Abstract 

Context. Close white dwarf - red dwarf binaries must have gone through a common-envelope phase during their evolution. DE CVn 
is a detached white dwarf - red dwarf binary with a relatively short (~8.7 hours) orbital period. Its brightness and the presence of 
eclipses makes this system ideal for a more detailed study. 

Aims. From a study of photometric and spectroscopic observations of DE CVn we derive the system parameters which we discuss in 
the frame work of common-envelope evolution. 

Methods. Photometric observations of the eclipses are used to determine an accurate ephemeris. From a model fit to an average low- 
resolution spectrum of DE CVn we constrain the temperature of the white dwarf and the spectral type of the red dwarf. The eclipse 
light curve is analysed and combined with the radial velocity curve of the red dwarf determined from time-resolved spectroscopy to 
derive constraints on the inclination and the masses of the components in the system. 

Results. The derived ephemeris is HJD rain = 2452784.5533(1) + 0.3641394(2) x E. The red dwarf in DE CVn has a spectral type 
of M3V and the white dwarf has an effective temperature of 8 000 K. The inclination of the system is 86^° and the mass and radius 
of the red dwarf are 0.41 ± 0.06 M Q and 0.37^oqq 7 Rq, respectively, and the mass and radius of the white dwarf are 0.5Uqo2 M b an d 
0.0136tg;g$f R s , respectively. 

Conclusions. We found that the white dwarf has a hydrogen-rich atmosphere (DA-type). Given that DE CVn has experienced a 
common-envelope phase, we can reconstruct its evolution and we find that the progenitor of the white dwarf was a relatively low- 
mass star (M< 1.6M Q ). The current age of this system is 3.3 - 7.3 x 10' years, while it will take longer than the Hubble time for 
DE CVn to evolve into a semi-detached system. 

Key words. Stars: individual: DE CVn - Binaries: eclipsing - Binaries: close - Stars: late-type - White dwarfs - Stars: fundamental 
parameters 



1. Introduction 

Large gaps remain in our knowledge of binary stellar evolution 
that affect our understanding of not only evolved compact bina- 
ries, but also of phenomena such as supernovae type la explo- 
sions, the rate of neutron star - neutron star mergers, and the 
number of gravitational wave sources in our Galaxy. The poorly 
understood physics of the common-envelope (CE) ph ase results 
in co nsiderable uncertainty on the binary evolution (Paczvhski 
Il976l) . During the evolution of a binary, the more massive star 
turns into a giant. When the initial orbi tal period is small enough 
(<10 years, iTaam & Sandquis3l2QO O) the envelope of the gi- 
ant will encompass the secondary star. The secondary and the 
core of the giant will spiral in towards each other in a common- 
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envelope. When the envelope is expelled a close binary consist- 
ing of the core of the giant, which will evolve towards a white 
dwarf, and the un-evoly ed secondary star may emerge (see e.g. 
iNelemans & Tbu1l2005h . 

The co mmon-envelope phase is ex pected to be very short (< 
1000 years, Taam & Sandquist 2000) and is therefore virtually 
impossible to observe directly. To study the effects of this phase 
it is best to focus on objects that have most probably undergone 
a common-envelope phase in their past. These we identify with 
binary systems containing at least one stellar remnant where the 
current orbital separation is smaller than the radius of the giant 
progenitor (usually with orbital periods < 1 day). 

Eclipsing close binaries offer the greatest possibility of de- 
riving precise physical parameters of the stars. The masses, radii 
and orbital separations give insight into the binary evolution and 
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specifically tell us if a CE phase has happened sometime in their 
past. 

Some examples of detached, close white dwarf - low- 
mass ma in-sequence star (red dwarf) e clipsing binaries are: 
RR Cae dBruchll 19991). NN Ser jHaefneJI 19891). EC134 71-1258 
(lO'Donoghue et al.ll2003l). G K Vir dGreen et alJl!986l) and RX 
J2130.6+4710 dMaxted et alj|2004l) . For a review o n deta ched 
white dwarf - red dwarf bi naries see e.g. iMarshl d2000t) and 
ISchreiber & Gansicke (12003b. The latest list of white dwarf - red 
dwarf binaries is in iMorales-Rueda et al.l d2005l with ten new 
systems compared to iMarshf 2000). It is necessary to study as 
many of these systems as possible to be able to compare their 
characteristics with population synthesis models and to find their 
space densities as a function of composition (e.g. white dwarf 
temperature, spectral type, age). 

DE CVn (RX J 1326.9+4532) is a relatively unstudied, bright 
(V = 12.8) eclipsin g binary. It was f irst discovered as an X-ray 
source by ROSAT ( Voge set al.l 19991) and has a proper motion of 
-0.198 ±0.002" yr" 1 in right ascension and -0.178 +.0.003" yr" 1 
in de clination as given in the USNO-B1 catalog dMonet et al.l 
120031 

This object was first studied photometrically by 
Robb & Greimel (119971) . From the light curve and the un- 
equal maxima they derived an orbital period of 0.364 days. The 
asymmetry in their light c urve needed a sta r spot to accurately 
model the light curve. Robb & Greimel (1997) measured 
eclipse depths of 0.054 + 0.010 magnitude in the R band and 
0.1 28 + 0.029 magnit u de in t he V band. 

iHolmes & Samusl (1200 ll) obtained UBVRI photometry for 
five nights in June 2000. They confirm the dependence of the 
eclipse depth with colours and found minimum depths of the 
eclipse of 0.10 magnitude in /, 0.15 in R, 0.30 in V, 0.60 in B, 
and 1.00 in U. The differences with wavelength band indicate 
that the two stars have very different colours. 

We note a differ e nce i n the eclipse depths as quoted 
by IHolmes & Samusl (f2 001) compa red to the value s from 
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give their values as being eclipse depths, when looking at the 
light curve we suggest that they have taken the difference be- 
tween minimum and maximum light instead of the difference 
between the start and m inimum of the eclipse which is used by 
Rob b & Greimel (1 19971) and in the present work. 

DE CVn consists of an M-type star with a spectroscopically 
unseen companion, presumably a white dwarf. Throughout this 
paper we will refer to the M dwarf as the secondary component 
and the probable white dwarf as the primary component. 

In Sect. [2] we describe our observations and reductions. The 
results are shown in Sect. [3] and the conclusions are given in 
Sect. E] 



Table 1. Log of the photometric data of DE CVn. UVic is the 
automatic 0.5-meter telescope of the Climenhage Observatory 
in Victoria, Canada. DAO is the 1.8-meter telescope of the 
Dominion Astrophysical Observatory. MDM is the 1.3-meter 
telescope of the Michigan-Dartmouth-MIT Observatory in 
Arizona. ULTRACAM are observations with this instrument at 
the WHT. T is the integration time per observation in seconds 
and # is the number of observations. 



Date 


Tel 


Filter 


T 


# 


1? Anr 1997 

i _ n lii ±77/ 


UVic 


ft 


99 


173 


71 Anr 1997 

Z. 1 n LII ±77/ 


UVic 


ft 


99 


151 


?? Anr 1997 


UVic 


ft 


120 


179 


94 Anr 1Q97 


UVic 


ft 




l?^ 


? Mav 1997 

_ ivxtxy ± 77 / 


UVic 


ft 


99 


175 


8 Mav 1997 

O lVJ.tlY ± 7.7 / 


UVic 


ft 


99 


76 


8 Mav 1997 

vj ly±a y 177/ 


UVic 


y 


120 


76 


9 Mav 1997 

7 ivitiy 1 77 / 


UVic 




^4. 
jt 




10 Mav 1997 


UVic 


clear 


33 


484 


1 Jul 1998 


UVic 


R 


140 


101 


7 Mar 2000 


UVic 


clear 


99 


52 


21 May 2001 


UVic 


R 


120 


169 


26 May 2001 


DAO 


B 


40 


110 


2 Jun 2001 


UVic 


clear 


40 


130 


21 Jan 2002 


MDM 


B 


60 


193 


24 Jan 2002 


MDM 


BG3S 


4 


597 


30 May 2002 


DAO 


B 


30 


39 


1 Jul 2002 


MDM 


BG38 


5 


216 


3 Feb 2003 


MDM 


B 


10 


350 


22 May 2003 


ULTRACAM 


u'g'i' 


1.3 


-750 


24 May 2003 


ULTRACAM 


u'g'i' 


1.3 


-2860 


25 May 2003 


ULTRACAM 


u'g'i' 


1.3 


-3780 


19 Jun 2003 


MDM 


B 


10 


166 


4 May 2004 


ULTRACAM 


u'g'i' 


1.3 


-13000 


5 Apr 2006 


UVic 


clear 


33 


359 



To reduce dead-time for the MDM observations, the CCD 
was binned 2x2, resulting in a scale of 1 .02" per binned pixel, 
and only a subregion of 256 x 256 (binned) pixels was read. 
These observations were reduced with standard packages in the 
Image Reduction and Analysis Facility (IRAFfl We derived dif- 
ferential photometry with respect to one comparison star (RA = 
13:26:59.6, Dec = +45:33:05, J2000) that is bright enough in the 
sparse field. The UVic and DAO observations were reduced with 
standard packages in IRAF. 

The ULTRACAM data were reduced with standard aperture 
photometry. Differential photometry was obtained with respect 
to two comparison stars located at RA = 13:26:28.08, Dec = 
+45:33:11.6 (J2000) and RA = 13:26:39.2, Dec = +45:34:56.1 
(J2000). The coordinates of DE CVn are RA = 13:26:53.2, Dec 
= +45:32:46.1 (J2000). 



2. Observations and reductions 

2.1. Photometry 

Our photometric dataset consists of various observations taken 
on a number of telescopes. Table Q] lists an overview of our 
photometric datasets taken with the 1.3-meter telescope of the 
Michigan-Dartmouth-MIT (MDM) Observatory (Arizona), with 
the 4.2-meter Willi am Herschel Telescope ( WHT) on La Palma 
with ULTRACAM dDhillon & Marshll200ll) , with the automatic 
0.5-meter telescope of the Climenhage Observatory in Victoria, 
Canada (referred to as UVic) and with the 1.8-meter telescope 
of the Dominion Astrophysical Observatory (DAO) located in 
Victoria, Canada. 



2.2. Spectroscopy 

A log of our spectroscopic data set is given in Table [2] The TLS 
echelle spectra are reduced with standard packages in IRAF. 
We could not correct the TLS spectra for flat fielding, because 
the sky/twilight flat fields would introduce only more lines in 
the spectra. The wavelength calibration was done with -300 
lines in the Th - Ar arc spectra with a root-mean-square resid- 
ual of -0.0036 A. We checked the stability of the spectrograph 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Table 2. Log of spectroscopic observations where we have 
given the date, telescope (Tel), the wavelength range (A) 
in A, the integration time (T) in s and the resolution (R) 
in A as derived from the FWHM of the arc lines. TLS 
= 2m telescope of the Thiiringer Landessternwarte 'Karl 
Schwarzschild' Tautenburg Echelle Spectrograph, Modspec = 
2.4m Hiltner Telescope (MDM) Modular Spectrograph, CCDS 
= 2.4m Hiltner Telescope CCD Spectrograph, Mklll = 2.4m 
Hiltner Telescope Mark III Spectrograph, DAO = Cassegrain 
Spectrograph at the 1.8m of the DAO. 





Date 


Tel 


A 


T 


# 


R 




11 May 1998 


TLS 


4300-5300 


900 


21 


0.14 


12- 


18 May 1998 


TLS 


5650-10 000 


900 


75 


0.16 


19 


- 22 Jan 2002 


Modspec 


4250-7550 


120 


62 


3.5 




16 Feb 2002 


Modspec 


4250-7550 


120 


5 


3.5 




19 Feb 2002 


Modspec 


4250-7550 


120 


10 


3.5 


5 


- 7 May 2002 


CCDS 


4180-5100 


300 


18 


3.2 




13 May 2002 


CCDS 


4180-5100 


360 


1 


3.2 




13 May 2002 


CCDS 


4180-5100 


600 


1 


3.2 




12 Jun 2002 


Mklll 


4250-7550 


90 


1 


3.7 




13 Jan 2004 


Modspec 


4250-7550 


180 


1 


3.5 


17 


- 18 Jan 2004 


Modspec 


4250-7550 


180 


2 


3.5 




6 Mar 2004 


Modspec 


4250-7550 


180 


2 


3.5 




30 Jan 2006 


ISIS 


3000-8000 


180 


1 


5 




30 Jan 2006 


ISIS 


3000-8000 


240 


6 


5 




15 May 2006 


DAO 


3500-5150 


720 


6 


4.8 



by using sky lines and these are not shifted between the ob- 
servations. Because of the unavailability of sky background or 
flux standards the echelle spectra are not sky subtracted or flux- 
calibrated. These spectra were taken during grey time, so we do 
not expect problems with the background level. 

The MDM (using the Modspec, CCDS, Mklll spectro- 
graphs), DAO and ISIS spectra were all reduced with standard 
packages in IRAF and these are, except for the DAO spectra, 
approximately flux calibrated (a seeing matched slit width was 
used, limiting the accuracy of the photometric calibration). 

3. Results 

3.1. Eclipse light curves 

T he first ephemeris of the primary eclipse in DE CVn is given 
in Robb & Greimell(ll997l) . We take the mid-eclipse times as the 
mid-point between the start of the ingress and the end of the 
egress. From o ur photometric observations and the times of min- 
ima given in lTas et al] (|2004|) we have determined the ephemeris 
of the eclipse minima by fitting a strai ght line to the cycle num - 
bers as derived from the ephemeris of Robb & GreimeJ d 1997b : 

HJD min = 2452784.5533(1) + 0.3641394(2) x E (1) 

No significant aliases were found near this period. The uncer- 
tainties (last digits) are derived for Ay 2 = 1 when we scale 
the individual errors to obtain that the reduced \ 2 — 1 ■ The 
value of the orbital period is further confirmed by the radial 
velocity analysis described in Sect. 13.51 The availability of two 
eclipses separated by only two cycles (-6109 and -6107) leaves 
no cycle count ambiguity. Two mid-eclipse times of Tas et al. 
(2452411.3156 and 2452412.4078) were rejected due to the 
large phase shifts (A<p > 0.01, in contrast to an average scat- 
ter of 0.004 f or all other eclipses) wit h respect to the updated 
ephemeris of Ro bb & Greimell 0997) leaving 23 mid-eclipse 
times for determining the ephemeris. When calculating the phase 



Table 3. Times of mid-eclipse for the binary DE CVn are given 
together with the cycle number and the corresponding time dif- 
ference. The numbers in parentheses a re the uncertainties on th e 
mid-eclipse times in the last digits. ° Robb & Greimel (1997); 
b UVic photomet ry; DAO photometry; d MDM photometry; 
" iTas et al.1 d2004l) ; f ULTRACAM photometry; « Not used for 
determining the ephemeris. 



HJD-2450000 


cycle 


A (HID) 


550.9214(16)" 


-6134 


-0.00060 


560.0243(16)° 


-6109 


-0.00118 


560.7531(20") 


-6107 


-0.00066 


562.9374(22)" 


-6101 


-0.00120 


570.9497(14)" 


-6079 


0.00004 


576.7749(14)" 


-6063 


-0.00100 


578.9606(6)" 


-6057 


-0.00013 


995.9003(18)* 


-4912 


-0.00010 


1620.7628(15)* 


-3196 


-0.00089 


2050.8125(16)* 


-2015 


0.00012 


2055.9103(7)" 


-2001 


-0.00003 


2062.8289(7)* 


-1982 


-0.00008 


2295.8782(4/ 


-1342 


-0.00003 


2298.7914(l) rf 


-1334 


0.00006 


2411.3156(1)"* 


-1025 


0.00517 


2412.4078(22)"« 


-1022 


0.00495 


2413.4958(4)" 


-1019 


0.00053 


2424.7839(6)" 


-988 


0.00031 


2456.8279(4/ 


-900 


0.00004 


2673.85488(14/ 


-304 


-0.00009 


2705.5359(3/ 


-217 


0.00079 


2727.3837(4)" 


-157 


0.00023 


2784.553370(25/ 





0.00000 


2809.6788(1)'' 


69 


-0.00019 


3830.7256(4)* 


2873 


-0.00040 



T 


I 1 I 1 I 
i' 


1 1 






/ 




, V, 

u 


o 








,-v-t - 

1,1,1 


1 , 1 



I I I I I I I I I 

-0.04 -0.02 0.02 0.04 



Phase 

Figure 1. Primary eclipse observed with ULTRACAM. 



difference with our new period, it turns out that these two data 
points have indeed a large phase difference (A<p > 0.01), so most 
probably the published times of minima are incorrect. The times 
of minima together with the cycle number and the corresponding 
time difference are given in Table [3] 

We observed a primary eclipse simultaneously in u' , g' and 
i' with ULTRACAM. These data show the largest eclipse depth 
of 1.11 ± 0.04 magnitudes in u' . The eclipse depths in g' and 
i' are 0.235 ± 0.004 and 0.028 + 0.004 magnitudes respectively, 
where we have fitted straight lines to the out of eclipse points and 
in eclipse points to derive these values. The difference is taken 
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as the eclipse depth. The uncertainties are derived for Ax 2 = 1 
when we scale the individual errors to obtain that the reduced 

3.2. Apparent magnitudes 

In Table |4] we give the magnitudes of DE CVn in and out of 
eclipse of the ULTRACAM data as obtained with the second 
comparison star as given in Sect. [XT] Also the MDM magnitudes 
and the magnitudes of the comparison stars are given. 

The MDM photometry taken on January 21, 2001 was ob- 
tained near phase 0.7, well outside the eclipse. The transforma- 
tion to standard magnitudes was derived using observations of 
four Landolt standard fields. Except for U - B, the transforma- 
tions derived from the standard stars all had scatter below 0.03 
mag. For U - B the scatter is below 0.08 mag. 

We derived the u', g' and i' magnitudes of DE CVn in and out 
of eclipse in the ULTRACAM data of May 24, 2003. The part 
just before the actual eclipse is taken for the out of eclipse mag- 
nitudes. All the ULTRACAM magnitudes were derived by using 
the two compariso n stars in the field as mentioned in Sect. 12.11 
Smit h et al] (120021) values were used to derive the u', g' and i' 
magnitudes. The magnitudes for DE CVn were the same us- 
ing either comparison star within ~0.01 mag in g' and i' and 
~0.02 mag in u' (lcr uncertainty). 

3.3. The nature of the components 

DE CVn has not been studied spectroscopically before. An av- 
erage low resolution spectrum taken with the ISIS spectrograph 
on the WHT on La Palma is shown in Fig. [2] Clearly visible are 
the absorption bands of TiO indicating an M-type star. Emission 
lines of Ha, HB, Hy, H<5 and Call H&K are visible as well. 

DE CVn is a single-lined spectroscopic eclipsing binary. We 
do not see any spectral features of the white dwarf in the over- 
all spectrum. The six low dispersion DAO spectra referred to 
in Table [2] were observed consecutively before, during and after 
an eclipse of the white dwarf by the red dwarf. The sum of the 
two spectra taken during the eclipse were then subtracted from 
the sum of the two spectra taken immediately before the eclipse. 
The resultant smoothed spectrum is plotted in Fig. [2] Using the 
spectra taken after the eclipse resulted in a similar spectrum. The 
strong hydrogen absorption lines are typical of a DA white dwarf 
and by visually comparing our WD spectrum with the ones in 
Wesemae l et al.l d 19931) we come to a spectral type of DA7 + 0.5 
which corresponds to a temperature of 7500 + 1000 K. The lack 
of residual Ca II K and narrow hydrogen emission lines gives us 
confidence that the subtraction was done correctly and the spec- 
tra did not need to be scaled. 

To determine the characteristics of the two components we fit 
the averaged ISIS spectrum with a composite model consisting 
of a white dwarf and a red dwarf. We first corrected the spectra 
for the radial velocity variations as a function of phase before 
averaging the ISIS spectra. 

The comparison spectra that were used to fit the data con- 
sist of a white dwarf model spectrum with a hydrogen atmo- 
sphere and temperatures between 1 500 and 17 000 K (kindl y 
made available to us by R Bergeron: Berge ron et all 199 111 19951) . 
From the most likely white dwarf mass and radius as derived 
from the eclipse fitting in Sect. |3.8| we derive a surface gravity 
for the white dwarf of logg ~7.8. Therefore we have used only 
white dwarf template spectra with surface gravity log g — 7.5. A 
red dwarf template (MOV to M6 V in integer types) together with 



the corres ponding absolu te visual magnitude was taken from the 
library of Pickles (1998). We first scaled the individual spectra 
to 10 pc before adding them together. We calculated the reduced 
X 1 of the fit to the average ISIS DE CVn spectrum for all the 
different model composite spectra to determine the nature of the 
components in this binary. The composition of the model spec- 
trum with the lowest reduced x 2 is taken as the best combina- 
tion. A more extended description of the fitting procedure will 
be given in a future paper. 

For the fitting of DE CVn we have excluded the wave- 
length regions around the emission lines of He, HB, Hy, H6 
and Ca II H&K and the earth atmosphere bands. The best fit 
consists of the combination of a red dwarf with spectral type 
M3V and a white dwarf with a temperature of 8 000 K although 
the corresponding formal reduced^ 2 is high (511). The second 
best fit has a Ax 2 = 50. When we take all the combinations 
with a x 2 < 1 000 the spectral type of the red dwarf stays the 
same, while the temperature varies between 7 000 and 9 000 K. 
Therefore we take an uncertainty on the temperature of 1 000 K. 
If we used \ogg = 8.0 instead of log g = 7.5 the spectral type of 
the secondary stayed the same, while the temperature changed 
to 10000 + 1 500 K. Combined with the temperature estimate 
from the eclipse (end of this section) and the difference spec- 
trum (see above) we decided to use the results from fitting with 
the logg = 7.5 models. 

The spectrum of DE CVn together with the best fit is shown 
in Fig. [2] Consistent results are found when fitting the MDM 
spectra which cover a shorter wavelength range. The discrepancy 
between the data and the fit of the spectrum are likely due to flux 
calibration errors, different intrinsic properties of the red dwarf 
such as metallicity, and the non-removal of telluric features in 
our spectra. 

Another method t o derive the seco ndary spectral type is by 
using the Ti05 index. Re id et al.l (119951) have given the best-fit 
linear relation between the spectral type of a late-type main- 
sequence star and its Ti05 band strength: 

S p = -10.775 xTi05 + 8.2 (2) 

where Ti05 is the band strength as defined by F w /F cont with F w 
the flux in the 7126-7135 A region and F cont as the flux in the 
7042-7046 A region. 

Using this definition we have calculated the Ti05 band 
strength and the corresponding spectral type for all the MDM 
spectra covering the wavelength range given above. The white 
dwarf contribution to this part of the spectrum is very small 
(< 5%), so this will not affect the ratio. The corresponding aver- 
age S p is 2.17, corresponding to spe ctral type M2. 

From the single M3 spectrum of lPicklesI d 19981) we have cal- 
culated the Ti05 band strength and corresponding spectral type 
as well. This gives a strength of 0.55 and a spectral type of 2.27. 
By comparing the value of t he M3 spec trum with the intrinsic 
variation as seen in Fig. 2 of Rei d et all d 1995b we see that this 
value is consistent with the values derived for DE CVn. From 
this we see that the different method s are consisten t with a red 
dwarf of spectral type M3 as used bv lPicklesI d!998l) . 

From the apparent magnitudes of DE CVn in eclipse and out- 
side eclipse we can derive the colour of the unseen white dwarf 
that is being eclipsed. This gives (u' - g')wD = 0.52 + 0.01 and 
ig' - i')wo = -0.26 + 0.02. By comparing this with the val- 
ues from the white dwarf models the u' - g' colour indicates a 
temperature of 7 000-9 000 K, while the g' - i' colour indicates 
6 000-8 000 K. This is fully consistent with the spectral mod- 
elling. 
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Table 4. u', g' and i' magnitudes of DE CVn and the two comparison stars as derived from the ULTRACAM data taken on May 24, 
2003. The uncertainties in u', g' and i' are 0.01, 0.01 and 0.02, respectively. The magnitudes and colours as derived from the MDM 
photometry are given as well with the uncertainty on the last digits in parentheses. 





u' 


s' 


j' 
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U-B 


B — V 


V-I 


DE CVn in eclipse 


16.43 


13.74 


11.65 










DE CVn out-of-eclipse 


15.31 


13.50 
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0.070(15) 


1.263(5) 


2.244(2) 


comparison 1 


15.16 
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comparison 2 
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Figure 2. A combined ISIS spectrum of DE CVn (black line) together with a composite template of an M3 V star from iPicklesI 
(1998) and a DA white dwarf with a temperature of 8 000 K and log g — 7.5 (grey/green line). Emission lines of Ha, Hfi , Hy and 
H6 are visible in the spectrum as well as the Ca II H&K emission lines and the TiO absorption bands of the M dwarf. No spectral 
line signatures of the primary are visible. The difference between the in and out of eclipse spectra shows the underlying white dwarf 
which is plotted in the upper left corner. 



3.4. Spectral line variations 



The spectra show emission lines of hydrogen up to H 10 and 
Call H&K emission. We have searched the normalized TLS 
echelle spectra for spectral lines showing radial velocity vari- 
ations, either in phase with the Balmer and Call H&K lines or 
in anti-phase. All lines identified are listed in Table [5] together 
with their equivalent widths (EW). No lines were seen to move 
in anti-phase with the Balmer lines. The lines that do not have 
an EW value are blended with sky lines so that we can not de- 
rive an accurate value for the EW. The EW of the bluest H lines 
were measured in the average ISIS spectrum for which we first 
removed the phase shifts in the individual spectra. 



3.5. Radial velocity curve 

Radial velocities of the Ha lines in the TLS spectra were deter- 
mined by fitting a single Gaussian line profile and a first order 
polynomial to the emission line and the surrounding continuum. 
The radial velocities of the H/3 and Hy lines in the TLS spectra 
were also measured in this way. The typical uncertainties on the 
radial velocities of the TLS spectra for Ha, H/3 and Hy are ~0.3 
km s _1 , ~0.4 km s _1 and ~0.7 km s _1 , respectively. 

The MDM spectra were cross-correlated with an M dwarf 
spectrum over the 6000-6500 A range. The uncertainties for 
these cross-correlated radial velocities are ~2.2 km s _1 . 

To derive the semi-amplitude of the radial velocity variations 
of the secondary and the systemic velocity we use the measure- 
ments of the Ha line in the TLS spectra, because these are the 
most accurate measurements. There is no spectral line feature of 



6 



E.J.M. van den Besselaar et al.: DE CVn: A bright, eclipsing red dwarf - white dwarf binary 



Table 5. Identified lines in the TLS echelle spectra and the ISIS (°) spectrum together with their equivalent widths. The lines that do 
not have an equivalent width measurement are blended with sky lines. A b s is the wavelength of the line as measured in the echelle 
spectra, while A is the wavelength corresponding to the given element which we identify with this line. The typical uncertainty on 
the EW is -0.1 A. 



/lobs (A) 


A (A) 


Element 


EW(A) 


/Us (A) 


A (A) 


Element 


EW(A) 


3797.20" 


3797.91 


H10 


-0.9 


6764.56 


6764.13 


Fel 


0.3 


3832.45" 


3835.397 


H9 


-1.8 


6768.72 


6768.65 


Til 


0.4 


3887.33" 


3889.05 


H8 


-2.1 


6773.05 


6772.36 


Nil 


0.3 
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4340.465 


Hy 
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6921.45 


6920.16 


Fel 


0.3 


4860.63" 
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H/3 


-4.1 
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7054 


TiO 


1.5 


5094.51 
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Fell 


0.4 
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0.6 
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Fel 
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0.3 
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1.3 
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blended 
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Fel 
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0.6 


6760.15 


6760.61 


Fel 


0.4 












Phase Phase 

Figure 3. Radial velocity measurements for the Ha lines. The Fi gure 4. Residuals of the Ho-, H# Hy and the MDM absorption 
solid line is the best fit to these velocities, as discussed in lines as compared with the best fit sinusoid to the Ha data. 
Sect. 13.51 The uncertainties on the points are smaller than the 
symbols. 



the primary star visible in the spectra, so therefore we cannot 
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Phase 

Figure 5. The equivalent width of the Ho- line in the TLS spectra. 
The variations are folded on the orbital period of DE CVn. A 
typical error bar is shown in the bottom left corner. 



derive the semi-amplitude of the white dwarf. We have used the 
function 



/(0) = r + ^2sin(27r(0-0 o )) 



(3) 



where K2 is the semi-amplitude of the secondary star and y the 
systemic velocity in km s to fit the data. The best fit gives 
y = -7.5 + 3 km s -1 , K 2 = 166 + 4 km s" 1 and <f> = -0.004 + 
0.004 for fitting only the TLS Ha variations. The uncertainties 
are derived from A% 2 = 1 when the reduced x 2 — 1 ■ The phase 
of the radial velocity curve is consistent with the time of mid- 
eclipse. The fit is shown in Fig. [3] together with the measured 
radial velocity variations of Ha. The residuals of Ha, H/?, Hy and 
the TiO absorption features are shown Fig. [4] When fitting for 
an eccentric orbit we find an insignificant eccentricity of e — 
0.02 + 0.02 so assuming a circular orbit is justified. 

No differences in the radial velocity curves of H/3, Hy and 
absorption with respect to the Ha observations are observed. 

3.6. Irradiation 

In close binary systems we may expect to see irradiation on the 
red dwarf due to heating of the close white dwarf. To investigate 
this effect we have measured the variations of the EW of the 
Ha line. In Fig. [5] we show the variations of the EW of this line 
folded on the orbital period. A typical error bar is shown in the 
bottom left corner of the figure. There is considerable variation 
of the Ha line, but this does not coincide with the time scale of 
the orbital period. During the eclipse the EW is larger compared 
to outside eclipse. If the Ha emission originates in the atmo- 
sphere of the secondary due to irradiation of the white dwarf, we 
would expect to see variations in the strength corresponding to 
the time scale of the orbital period. 

In the similar system RR Cae, consisting of a white dwarf 
of 7 000 K and an M6 or later sec ondary, a sim ilar effect of a 
larger EW during eclipse was seen (lBruchT l999). This was con- 
tribut ed to the emi ssion lines being intrinsic to the secondary 
itself (lBruchlll999l) . In analogy, we therefore conclude that the 
Ha emission in DE CVn is also due to activity in the chromo- 
sphere of the M dwarf in contrast to being caused by irradiation. 
Furthermore, we can rule out Ha emission emerging from an ac- 
cretion disk, because of its regular radial velocity profile and the 
single peaked lines. 




Phase 

Figure 6. Ti05 band strength of the MDM spectra as a function 
of orbital phase. The uncertainty on each point is smaller than 
the symbol. 




20 30 
Time (min) 

Figure 7. A flare of DE CVn observed with ULTRACAM. 



If the secondary star is irradiated by the white dwarf we 
would also expect a variation of the Ti05 band strength (which 
is discussed in Sect. 13.3b with respect to the orbital period due 
to heating in the atmosphere of the secondary. We see variations 
(see Fig.|6j, but not with the orbital period. Thus also in this case, 
we conclude that the variations are intrinsic to the secondary. 

3.7. Flare 

During one of the ULTRACAM observing runs (25 May 2003) a 
flare was observed in all three bands starting at about 23:23 UT. 
This part of the light curve is shown in Fig. [7] The observed part 
of the flare lasted ~39 minutes. The flare has a complex struc- 
ture with several peaks close together and a rebrightening during 
the decay part after the first peaks. Unfortunately the observation 
ended before DE CVn had returned to the quiescent (pre-flare) 
luminosity (most clearly seen in the u' band). From the observed 
count rates during quiescence and the average count rate dur- 
ing the observed part of the flare we calculated that DE CVn 
increased in brightness on average by 20%, 5% and 1% of the 
quiescent flux (so white dwarf + red dwarf) during this flaring 
period in u' , g' and 2' respectively. 

To model the flare, we have taken a blackbody spectrum 
added to the best fitted template spectrum (see Sect. 13.31 ). A 
blackbody might not be a good approximation, but it is used in 
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Figure 8. Temperature variation of the flaring area during the 
flaring period. The flaring area has a radius of 0.038/? Q . 



lower inclinations we find that the red dwarf either fills its Roche 
lobe or is larger than its Roche lobe. As there is no evidence for 
mass transfer in the system, the red dwarf must be smaller than 
its Roche lobe. The reduced^- 2 for these fits to the eclipse light 
curve are not significantly different, illustrating the fact that an- 
other constraint on the system is required for a unique solution. 

The semi-amplitude of the radial velocity variations of the 
secondary is derived from the spectra in Sect. 13.51 The spectral 
type of the secondary is determined in Sect. |3.3| by fitting model 
spectra to the spectrum of DE CVn. These parameters are taken 
as known input parameters in our two-circle model. 

From the spectral type of the secondary and assuming 
zero-age main-sequence masses and temperatures the mass of 
the red dwarf should be between 0.3 and O.5M (see e.g. 
Kirkpatri ck et al.l Tl 99 lb . For every possible inclination we in- 
put the mass of the secondary in steps of O.O1M between these 
values. Furthermore, we use Kepler's third law and the mass- 
radius relation for white dw arfs from Eggleton as quoted by 
IVerbunt & Rappaporj(ll988l) : 



other analyses of flare stars as well (de Jag er et ailll986l [l989). 
We have assumed a constant (but not pre-determined) flaring 
area that is facing us with only temperature variations over time. 
We have assumed temperatures between and 20 000 K in steps 
of 200 K and a constant radius for the flaring area. 

For every point in time we have taken the three correspond- 
ing observed magnitudes. We have derived the u', g' and i' mag- 
nitudes for our model by convolving the filter curves with our 
model spectrum for each temperature. For each point in time we 
have calculated the reduced^ 2 for each temperature and we have 
taken the temperature with the lowest reduced x 1 for this point. 
After this the average reduced y 1 ° ver the time points was 
calculated. 

We performed this method for several constant flaring areas 
with radii between 0.010fl o and O.O53/? with steps of 0.001fl o . 
The result with the lowest average reduced x 2 h as a radius 
O.O387? and the temperature variations are shown in Fig. [8] The 
area that flares corresponds to 0.9% of the visible stellar disk of 
the red dwarf. 

3.8. System parameters 

Because spectral features from the white dwarf are not visible 
we can only measure the radial velocity curve of the secondary 
star. The inclination of the system is constrained by the eclipses 
in the photometric light curve. 

We have fitted the ULTRACAM g' light curve with a simple 
model light curve that results from using two overlapping cir- 
cles representing a white dwarf and a red dwarf that is (partly) 
obscuring the white dwarf (see Van Ham et al. 2007, submitted). 
It is assumed that the light intensity is proportional to the visible 
part of the white dwarf. We have used the light curve in counts 
scaled to range between 1 (out of eclipse) and (in eclipse). This 
way we only derive the radii of the two stars with respect to the 
orbital separation (a, i.e. Rwd/q and Rrd/o). 

Limb darkening of the white dwarf could play a role. When 
we include limb darkening in our model the white dwarf radius 
increased by a few per cent. This is negligible with respect to the 
uncertainty in the inclination so therefore we have neglected the 
effect of limb darkening in our method. 

From our model we obtain the combinations of the white 
dwarf and red dwarf radii with respect to the orbital separation 
for inclinations between 75° and 90° with a step size of 1°. At 



fl WD = 0.0114- 



/Mwe>\-3 _ 
V M CH / 

M WD 



Men 
M WD 



V Mp J \ M P J 



(4) 



where 7?wd is in solar radii, Mwd in solar masses, Men 
= \AAM Q and M P is a constant whose numerical value is 
0.00057M o . 

For every inclination we obtain a value of Swd/i from the 
eclipse fitting. By combining Kepler's third law and eq. (01 we 
derive the mass and radius of the white dwarf for every inclina- 
tion and secondary mass. With the radius of the white dwarf and 
the result of the eclipse fitting we derive the orbital separation 
and the radius of the secondary star. 

By assuming a circular orbit, which is justified by the ra- 
dial velocity curve, we can constrain the inclination by using the 
mass function given in eq. © to derive the semi-amplitude of 
the red dwarf for every inclination and secondary mass, and the 
system parameters derived from the procedure described above. 
With these parameters and assuming a circular orbit we can use 
the following function for the mass function: 



M? sin 3 i , 
1 =(1.0361 x 10" 



(M 2 +Mi) z 



kk p 



(5) 



with Mi 2 in M Q , K 2 in km s _1 and P in days. 

We take a 3<x uncertainty on the measured semi-amplitude 
of the radial velocity variation of the secondary and compared 
these with the semi-amplitudes for each possible solution. We 
selected only those solutions that satisfied the radial velocity cri- 
terion. This constrains the inclination of DE CVn to > 82°. The 
combined radial velocity and eclipse constraints lead to system 
parameters where the mass of the white dwarf is between 0.48 
and 0.65M o and radius between 0.01 17 and 0.0140fi o . The red 
dwarf has a mass between 0.30 and O.5OM and radius between 
0.36 and O.51/? . The orbital separation of the system is between 
2.00 and 2.25fl . 

By comparing our white dwarf radii with the correspond- 
ing masses for car bon core or oxygen core white dwarfs from 
Pan ei et all (S000) our minimum and maximum mass for the 
white dwarf would be at most 5% larger. Furthermore, all the 
red dwarf masses that we used in deriving the system parame- 
ters are possible solutions. 
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In our ULTRACAM data set, we see hints for small out-of- 
eclipse light curve variations and variations between the differ- 
ent observations, but this data set does not cover the complete 
orbital period. In our eclipse fitting procedure we only use the 
section between phases -0.035 and 0.04, so here the effect of 
these variations can be neglected. 

3.9. Space velocity and distance 

The fitting procedure described in Sect. 13.31 also determines a 
distance to DE CVn of 28 + 1 pc. The proper motion is known as 
well (see Sect. Q} and the systemic velocit y is derived from fit- 
ting th e radial velocity curve in Sect. 13.51 Jo hnson & Soderbloml 
d 1987) give the equations to calculate the space velocities (U, 
V, W) and we use (U@, V@, W@) = (10.00 + 0.36, 5.25 + 0.62, 
7.17 + 0.38) in km s _1 as the velocity of the Sun with re s pect t o 
the local standard of rest as given bv lDehnen & Binnevl (11998). 
The space velocities of DE CVn, defined as 

(u, v, w) = (U, V, W) - (U Q , V q , W ) (6) 

are (u, v, w) = (-15.6 + 0.6, -40.7 + 1.5, -3.2 + 2.9) km s _1 . 
U is defined as positive in the direction of the galactic center, 
V is positive in the direction of the galactic rotation and W is 
positive in the direction of the galactic North Pole. The derived 
space v elocities for DE C Vn are consistent with being a thin disk 
object dFuhrmannl2004l) . 

4. Discussion and conclusions 

4.1. Binary parameters 

We have analysed various spectroscopic and photometric obser- 
vations of the eclipsing binary DE CVn. The light curves show 
total eclipses of the primary by the secondary star. 

The low resolution average ISIS spectrum is fit with model 
composite spectra to determine the spectral type of the M dwarf 
and the temperature of the white dwarf. The best fit is derived for 
an M3V secondary star and the temperature of the white dwarf 
(with log g = 7.5) of 8 000 + 1 000 K. This is consistent with the 
temperature estimates derived from the eclipse depths (Sect. l3.2l l 
and difference spectrum (Sect. 13.31 ). and the spectral type from 
the Ti05-index (Sect. [331 . 

From the echelle spectra taken with the TLS we have ob- 
tained radial velocity curves of Ha, Hfi and Hy and also TiO 
absorption bands from the MDM spectra. These curves show a 
semi-amplitude for the secondary of 166 ± 4 km s _1 . 

Combining the eclipse constraints with the radial velocity 
amplitude results in a set of solutions satisfying the constraints as 
described in Sect. 13. 81 Determining the most likely value with an 
uncertainty on it is very complex. Due to the interdependency of 
the relations used, calculating the uncertainties is not straightfor- 
ward. As an example, we show the histogram of the possible so- 
lutions for the white dwarf mass in Fig.|9]which shows an asym- 
metrical distribution. The average or the mean does not coincide 
with the peak of the distribution. We take the most likely value 
as the highest point in the distribution. To determine the uncer- 
tainty on the most likely value we have taken the value for which 
68% of the solutions are enclosed around the most likely value. 
The difference between these values is taken as our uncertainty. 
We have done this separately for each side of the distribution. 
For the white dwarf mass and radius this resulted in 0.5 1 ^ 02 ^0 
and 0.0136^oooo| R . In the same way, the distributions of the red 
dwarf mass and radius result in 0.41 + 0.06 M Q and 0.37^ ^ R Q 
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Figure 9. Distribution of the possible solutions for the mass of 
the white dwarf. The thick solid line (green/grey) is the most 
likely value that we have used. The thin solid lines (green/grey) 
are the values between which 68% of the possible solutions are 
located. The thick dashed line (red/darker grey) is the average 
value with the standard deviation in thin dashed lines (red/darker 
grey). The median of the distribution is shown as the dotted line 
(blue/dark grey) 



respectively. The orbital separation is 2.07 ^ R q and the in- 
clination is 86^ °- These most likely values are consistent with 
each other within the uncertainties. 

4.2. Progenitor 

From the current small orbital separation we conclude that there 
was most probably a CE phase in the past in which the orbit has 
shrunken significantly. At the onset of the CE phase, the initially 
more massive star has evolved to a star with giant dimensions. 
We assume that the core of the giant star at the start of the CE 
phase had the same mass as the present day white dwarf. This 
can be used to reconstruct the evolution of DE CVn to find the 
possible progenitors. 

To do this we take single main-sequence stars with masses 
of 1 to 8M in steps of O.1M . We evol ve these stars with the 
method described in iHurlev et al.l (l2000h . When the core mass 
(M c ) of the stars has reached the mass of the white dwarf that 
we observe in DE CVn (0.5 1M ) the evolution is stopped. Then 
we check if the radius of the star at this point (/?,•) corresponds 
to the largest radius during the evolution up to this point. We 
assume that the star fills its Roche lobe in the giant phase. 

The most likely mass of the white dwarf of DE CVn falls 
in a mass range where many progenitor stars cannot fill their 
Roche lobe with this core mass. The reason is that on the first gi- 
ant branch their core grows to ~O.48M , when the helium flash 
happens. Then the star contracts while the core mass still in- 
creases. When the star expands again to ascend the asymptotic 
giant branch (AGB), the core mass has become larger than the 
most likely white dwarf mass in DE CVn. However, this is very 
sensitive to the core mass: for a core mass of 0.5 1M we find 
only one progenitor, while for O.52M we find six. We therefore 
use O.52M , as a compromise between staying close to the most 
likely mass yet allowing as many progenitors as are allowed by 
the inferred mass range. 

The possible progenitors are given in Table [6] The first three 
columns in this table are the initial mass of the main-sequence 
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Table 6. Parameters of different white dwarf progenitors at the time the core mass reaches 0.52M o . Mo is the initial mass of the 
main-sequence progenitor. M, and R, are the mass and radius of the star when it has reached the required core mass. This means that 
the first five are AGB stars, while the last one is a red giant. f evo i is the time between the birth of the binary and the start of the CE 
phase, a, and P, are the initial orbital separation and initial orbital period, cce is the fraction of the orbital period used for ejection 
of the CE, while A is a numerical value dependent on the structure of the star. The time from main-sequence binary to a CV is given 
as t tot . 



M (M ) 


M, (M Q ) 


Ri (Re) 


4voi (yr) 


a, (R e ) 


P,(d) 


Q-CE/l 


tmt (yr) 


1.051 


0.923 


231.165 


6.51- 10" 


513.052 


1166.530 


0.03 


2.61-10 1 " 


1.191 


1.085 


224.726 


4.95- 10 9 


483.386 


1007.611 


0.05 


2.45-10 10 


1.319 


1.230 


218.584 


3.87-10 9 


459.140 


890.428 


0.08 


2.34- 10 10 


1.439 


1.365 


212.881 


3.10-10 9 


438.689 


799.406 


0.11 


2.27- 10 10 


1.555 


1.494 


207.045 


2.54- 10 9 


419.770 


722.411 


0.14 


2.21-10 10 


3.400 


3.399 


40.222 


2.74- 10 8 


71.425 


35.849 


5.83 


1.99-10 10 



progenitor of the white dwarf (Mo) and the mass and radius of 
the star at the time the evolution is stopped (M, and /?,•), The 
evolution time of the star until this point is given as f e voi- The first 
five possibilities are stars that reach the required core mass while 
on the AGB, while the last one is peculiar: it reached the required 
core mass when it had a non-degenerate core. If the progenitor 
was such a star, the system would have come out of the common- 
envelope phase as a helium-burning star (most likely observable 
as a subdwarf B star) with a low-mass companion. Only after 
most helium was burned to carbon and oxygen would the star 
have turned into a white dwarf. 

There are no possible progenitors with initial masses be- 
tween 1.6-3.4M Q . For mass > 1.6M the core mass of the giant 
on the red giant branch (RGB) has not yet grown to the required 
mass. After the RGB phase the star contracts again, while the 
core mass is growing. When the star reaches the AGB phase, 
the core mass is already larger than the required mass. This re- 
sults in no possible progenitor star. At the moment of sufficient 
core mass, the star needs to be larger than any time previously 
in its evolution otherwise the star would have started Roche lobe 
overflow earlier in its evolution when the core mass was not yet 
massive enough. 

For stars > 3.4M the core of the giant at the start of the 
RGB phase is already larger than the required core mass, leav- 
ing no possible progenitor for the present day white dwarf. The 
progenitor with an initial mass of 3.4M is just on the edge of 
being a possible progenitor for the present day white dwarf in 
DE CVn. 



4.3. Common-envelope phase 

We assume that a possible progenitor fills its Roche lobe at the 
st art of the CE pha se. By using the formula for the Roche lobe 
of lEggletod(fl983h 



Pi 



0.49^ 



0.6^3 + ln(l + qi) 



(7) 



(with q = M//M2) we can calculate the orbital separation (a,) 
(and thus the orbital period P,) at the moment the evolution was 
stopped. These values are given in Table|6]as well. 

From the initial orbital period (a,, Table |6]l and the present 
day orbital period, we note that a large orbital shrinkage has 
taken place. Together with the rather extreme initial mass ra- 
tios, this implies that DE CVn went through a common-envelope 
phase during its evolution. 



We take the parameters of the giants derived above as the 
ones at the start of th e common-envelope phase. With the use of 
(IdeKool et al. l ll987l) 



GMi(Mi-M c ) (GM c Mi GMiM 2 

— = ffCE^ 



Ri 



2d 



f 



2a, 



(8) 



we calculate ctce^ where A is a numerical value dependent on 
the structure of the star and c?ce the fraction of th e orbital energ y 
that is used for ejection of the common envelope (Ide Koo]|l 992). 
for every possible progenitor. <xqe is usually < acE < E but 
for do uble white dwarfs it is found to be > 1 ( Nelem ans & Toutl 
2005). The qtce^ values for the possible progenitor stars are 
given in Table [6] as well. 

The most mas sive progenitor has q-ce^ of 5.83. 
iDewi & Taurisl d2000t) give A values as a function of the 
radius of the giant for stars with M > 3M Q . From their table 
we conclude that for our possible progenitor of mass 3.4M A 
is around 0.6-0.9 indicating czce Z 6 which is very unlikely. 
Therefore we rule out such a star as the possible progenitor. 

The other possible progenitors have acE^ < 1 indicat- 
ing a rather inefficient ejection of the CE (e.g. a large frac- 
tion of the binding energy is not used for ejection of the CE). 
These values are in agreement with th e values found for other 
white dwarf - red dwarf bi naries (e.g. iNelemans & Tout! [2005: 
iMorales-Rueda et al]|2005l) . and we conclude that the progeni- 
tor star of the present day white dwarf was a main-sequence star 
with M < 1 .6M . 

4.4. Time scales 

Now that we can reconstruct the evolution of DE CVn we can 
determine the corresponding time scales. From the possible pro- 
genitors of the white dwarf we find that the time it took before 
the CE phase started is 2.5-6.5xl0 9 years. 

The temperature of th e white dwar f is 8 000 + 1 000 K. 
From the cooling tracks of Wood! ( 1 19951) as shown in Fig. 4 of 
ISchreiber & Gansickd (12003b we can determine a cooling age 
(?cooi) of the white dwarf of ~8 x 10 s years. The current age of 
the system is therefore 3.3-7.3xl0 9 years. 

Due to angular momentum loss, the system will evolve to- 
wards a semi-detached cataclysmic variable (CV) phase. To de- 
termine the time that it will take for DE CVn to become a semi- 
detached binary we first need to know the orbital period at which 
mass transfer starts (P s d). This period follows from the Roche 
geometry and Kepler's third law: 



Psd = 2tt 



GM 2 (i + f|)(^y 



0.5 



(9) 
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For DE CVn we derive an orbital period at the start of mass 
transfer of 3 hours. 

The orbital period will shrink due to loss of orbital an- 
gular momentum. For low-mass stars it is often assumed 
that the dominant mechanism is (d isrupted) magnetic breaking 
dVerbunt & Zwaan|[l98Tt lKindll988h . In short, disrupted mag- 
netic breaking occurs in a close binary when the tides force the 
secondary star to be co-rotating with the binary, while magnetic 
fields in the secondary force the stellar wind to co-rotate with 
the secondary star. When this exert a spin-down torque on the 
secondary star this must extract angular momentum from the bi- 
nary orbit. At a mass of ~O.3M the secondary becomes fully 
convective and therefore loses its dynamo (or at least changes it) 
so that magnetic breaking is no longer the dominant source of 
angular momentum loss. 

By assuming that the angular momentum loss is due to dis- 
rupted magnetic braking only until the secondary becomes fully 
convective at a secondary mass of at least 0.3M o , the time it wi ll 
take before mass transfer starts is (Schreiber & Giinsicke 2003): 



t 2 - 

hi ~ 



2.63- 10 29 Gm 



1 D-2 



(2n)-{M x +M 2 )- 



(10) 



where M\ and M2 are in M and R2 in R e . This corresponds to 
a time of 1.9 x 10 10 years before DE CVn becomes a CV, which 
is just longer than the Hubble time. 

Systems such as DE CVn will not contribute to the current 
sample of CVs, unless the loss of angular momentum in the cur- 
rent detached white dwarf - red dwarf phase is much higher than 
that g iven by magnetic braking alone (see e.g. iBrinkworth etaU 
2006, where an angular momentum loss mechanism ~100 times 
greater in strength than the currently accepted value seems to be 
required to explain the rate of period decrease in NN Ser). 
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